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1 Introduction 



The influence that boundary conditions have on different spectral functions 
is an active field of research. In quantum field theory, spectral functions 
of particular interest are the zeta function and the heat kernel. Their de- 
pendence on the boundary condition is well understood for a large variety 
of boundary conditions (for a recent review see pQ). However, for some 
boundary conditions an understanding of elementary properties of spectral 
functions is still lacking. This is the case for generalized (or chiral) local bag 
boundary conditions [21 El • 

These boundary conditions involve an angle 9, which is a substitute for 
introducing small quark masses to drive the breaking of chiral symmetry 
[SI lU E] • The influence the parameter 9 has on various correlators was ana- 
lyzed in detail in [Hj for the two dimensional Euclidean ball. In reference [£] , 
the heat kernel and the eta function were analyzed in the two dimensional 
cylinder. The situation of an arbitrary dimensional ball was considered in 
[7]. The results found in the above articles suggest that general properties 
of spectral functions, like the pole structure of the associated zeta func- 
tion and the form of the asymptotic expansion of the heat kernel, are the 
properties resulting from strongly elliptic local boundary conditions. For 
the generalized bag boundary conditions considered here, this property has 
not been proven. In fact, strong ellipticity is not clear at all because, for 
9/0, the boundary conditions are of mixed oblique type and, under 
certain circumstances, oblique boundary conditions are not strongly elliptic 
[HI El- However, after introducing some basic notation and properties of the 
boundary conditions, we will prove in Sections 3 and 4 that generalized bag 
boundary conditions are indeed strongly elliptic boundary conditions. Based 
on this observation, in future investigations one might envisage a determi- 
nation of heat kernel coefficients for these boundary conditions. In recent 
years, a conglomerate of methods has been proven to be very effective in 
the determination of this asymptotics 1 . Among the methods are special 
case considerations, which form the basis of the second half of our paper. 
In particular, we will determine the local heat kernel and the zeta function 
on cylindrical product manifolds. Results are given in terms of boundary 
data, much in the way it is possible for spectral boundary conditions |1U| . 
The Conclusions will summarize the main results and describe their possible 
future applications. We refer to [11-20] for other important related work in 
the field. 
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2 Basic properties of chiral bag boundary condi- 
tions 



In this section, we will establish some notations for the problem at hand, 
i.e., the Euclidean Dirac operator acting on spinors satisfying local (chiral 
bag) boundary conditions 

Let m = 2rh be even and let P = «7jVj be an operator of Dirac type on 
a compact oriented Riemannian manifold of dimension m. Let V denote the 
spinor space; dim(V) = 2 m . An explicit representation of the 7-matrices 
is provided in Appendix A. They are self-adjoint and satisfy the Clifford 
anti-commutation relation ()A.1|) . Near the boundary, let e m be the inward 
unit normal and let 7 m be the projection of the 7-matrix on the inward 
unit normal. In addition let 7 be the generalization of '75' to arbitrary even 
dimension, 7 = (— i) m 7i...7 m . 

We set 

X = iie 8l lm 
and use the relation 7^7 + 77^ = to compute 

2 -07 - 07 ~ 07- -07 -, 

X = -je 7 7m76 Ira = 7 7e 7 7 m 7 m = 1. 

We define 

n± :=\{l± X ) 

and have 

n| = n± and n_n+ = n+n_ = o. 

Note that these two projectors are not self-adjoint (except for the par- 
ticular case 9 = 0). Rather, calling their respective adjoints IT* and IPL, 
one has ILj- := |(1 ± ije~ e ^j m ), and the following equations hold 



n* n+ = cosh (6*7) exp (-6» 7 )n H 
n* n + = sinh (07) exp (— 6»7)LT 4 
n* n_ = cosh (6*7) exp (-6»7)n_ 
n^n_ = sinh (#7) exp (-07)n_ 



: II* cosh (07) exp (—07) 
II* sinh (07) exp (—07) 

: II* cosh (07) exp (—07) 
11*^ sinh (07) exp (—07) . 



We use II_ to define boundary conditions for P. Similarly, we shall let 
B := II_ © U_P define the associated boundary condition for P 2 . 

In the following two sections, we will show that (P, Il_) and (P 2 ,B) 
define strongly elliptic boundary conditions and, as a result, we can assume 
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standard results on the meromorphic structure of eta and zeta-invariants 
hold |21| . Otherwise stated, we will prove statements (1) and (2) of the 
following Theorem 

Theorem 2.1 

(1) (P, Il_ ) is strongly elliptic with respect to C — R + — R . 

(2) (P 2 ,£>) is strongly elliptic with respect to C — R + . 

(3) (P,n_) is self-adjoint. 

(4) (P 2 ,B) is self-adjoint. 

Statements (3) and (4) are well known to hold 3 and so we concentrate on 
statements (1) and (2). 



3 Ellipticity of the first order boundary value prob- 
lem 

Proof of (1): We use Lemma 1.11.2 (a) of (2^ to prove assertion (1). Note 
that the special case 9 = defines standard mixed boundary conditions 
and the theorem is known to hold for this case. Let x = (y\, y m -i, x m ) 
be coordinates near the boundary where x m is the geodesic distance to 
the boundary and where y = (y\, y m -i) are coordinates on dM. Let 
Udy a £T*{dM). 

Following [2TJ we define 

A) = -* 7m I lata - A J for (£, A) / (0, 0) and A R - {0}. 

\a<m / 

We then have (?(£, A) 2 = (|£| 2 - A 2 )l. As (|£| 2 - A 2 ) iR, we may let Vj be 
the span of the eigenvectors of q(£, A) with positive/negative real parts. We 
let 

W := Kernel(n_) = Range(n+). 

Using Lemma 1.11.2 (a) of |2j, we prove assertion (1) by verifying that LI_ 
is an isomorphism from V^(£, A) to W = Range(n_). This is equivalent to 
showing 

vl(n,\)nw = {0}. (3.1) 

We change variables slightly setting A = — i/j, where /i iH — {0} and 
^ (0,0). We then have 

M) = "'7m ^ lata + M7m and /l) 2 = (|£| 2 + 

a<m 
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The next step in the proof is to reduce the problem to a collection of effective 
two-dimensional ones. We make use of the properties of the 7-matrices. First 
note that the elements 

Tl := 27273, ...,Tm-l := «7m-27m-l, 

mutually commute and, in addition, they commute with 71 and j m . Thus, 
for j, k = 1, rh — 1, we have 

TjTk = TkTj, Tj = 1, 

Tjjl = JlTj, Tj<y m = J m Tj. 

So we can choose a set of simultaneous eigenvectors of Tj with eigenvalues 
Pj = ±1. We denote by p = (pi, Pm-i) the collection of simultaneous 
eigenvalues of Tj and we define the associated simultaneous eigenspaces by 

Vp = {v £ V : nv = piv}. 

The vector space Vp is preserved by 71, 7 m and 7. We use this fact to 
decompose Vp into its chiral parts, 

Vp = Vf Vp where = {v G Vp : 71; = ±v}. 

Let yJ" = span-ft;^}. Then since 7 m 7 = —77™ we have 

span{v p,-f m vp}. 

The vector spaces provide the decomposition 

V = @Vp, dim^=2, 

p 

and the problem completely decouples into two-dimensional spaces. 

On Vp one easily computes, using the definition e(p) = p\ X ... x pm-i, 

7l = e( ^(-°i 0)' 7m= (i 0)' 7= (o -1)' 

thus reproducing the two-dimensional Pauli-matrices up to the standard sign 
ambiguity. 

We note that the kernel of ITL is determined by the eigenvectors of x- 
Thus to establish (|3.1|) we shall need explicit representations of x an d q(£, p) 
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acting on Vp. To calculate q(£, /i) it is possible to choose coordinates such 
that £2 = ■•• = £m-i = 0. It is then easy to see that 

The eigenvectors of x follow, 

Q \ I Q 



-le " / \—ie 



and we compute 



le J \fj,g + ie(p)£,\e y 

where g = ±1. Assertion (1) holds if these are not multiples of each other, 
that is if 

det ( ~ f e Zl ) = -2ie(p)g^ + e"") + 0. 

V fJ,g + ie[p)i;ie —le J V / 

Since the first term is purely imaginary and /i zR— {0} this proves assertion 
(•)• 

In fact, to prove what we announced (i.e., that Il_ is an isomorphism 
from V"_5(£,A) to W = Range(II_)), it is enough to consider only g = +1. 
But since VH{— £, —A) = Vf.(£, A), considering p = ±1 gives the same condi- 
tion. 

Although assertion (2) follows from assertion (1) and Lemma 1.11.2 (b) 
in |21j . we prefer to give a second proof showing that the boundary opera- 
tor B involves tangential derivatives. This makes apparent that the chiral 
boundary conditions are non-standard boundary conditions. 



4 Ellipticity of the second order boundary prob- 
lem 

When considering spectral properties of the square of the operator of Dirac 
type, P 2 = («7jVj) 2 , the boundary condition imposed through B is 

n_v| 9A f = 0, (4.i) 

n_7 i V^| aA f = 0. (4.2) 
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The second boundary condition Q4.2J1 can be rewritten as an oblique bound- 
ary condition involving tangential derivatives. To this end note 



with the tangential 7-matrices j a , where a = 1, 
boundary condition (|4.2jl to be written as 



, m 



1. This allows the 



= n_(-i7jVj)V' \dM = «7mn_(— iTjVjOV'laAf 

= 7 m II_(7 m V m + 7 a V a )V> lajvf = (IT* V m + 7 m 7 a II* V a )V> \dM 

= (n;v m + 7m7a n*v a )(n_ + n + )^| 9M , (4.3) 

where V m is the interior normal derivative. The boundary condition con- 
tains tangential derivatives and the conditions imposed through B could 
thus be termed of mixed oblique type. 

Proof of (2): To study the ellipticity of the boundary value problem, 
we introduce the "partial" leading symbol of P 2 



&L(y,Xm,u,-id m ,X) 
and the graded symbol, a g , of B 



<7„ 



-7 a w a n* 



o 



Strong ellipticity requires that the problem 

—id m , A)*(y, x m ,u, A) 



with 



and 



(To 



r=0 



n_a 

Z7 m ni9 m a 



(4.4) 
(4.5) 

(4.6) 



r=0 



has an unique solution. 

Now, the solutions to (|4,4[) and 1)4. 5 Jl are 



W(x m ,u,\) = * exp(-Ax m ) 



7 



where A = +Vuj 2 - A. Note that »(A) > for A e C - R + . 
The condition (|4.6|) . when applied to them, reads 

/ n_ o \ / * \ = ( n_a \ 
V - 7a o; n* i 7m n* y V -a* y V -i 7m n*Aa ) ■ 

This gives a system of two equations. After multiplying the second one by 
i 7m , and using = n_\I>o + n + ^o 5 one obtains 

n_* = n_a, (4.7) 

and 

{-i~im{laOJaW_Tl + + All* 11+) tf = 

An;« - (-i 7m ( 7aWa )rrin_ + An + n_) ^ • (4.8) 

We use ((TTj) and substitute (|4~7)) into (|4~%)) to see 

exp (-0 7 ) [-i^mila^a) sinh (0 7 ) + A cosh (0 7 )] 11+ # = 

All* a - [-i 7m ( 7a w a ) cosh (0 7 ) + A sinh (0 7 )] exp (-0 7 )II_a . 
So, the problem has an unique solution if the matrix 

M = -i 7m ( 7a Wa) sinh (0 7 ) + A cosh (0 7 ) = A sinh (0) + A cosh (0 7 ) 

is nonsingular, where we introduced A = — i 77m7a w . 

But A* = A; so, it is diagonalizable. Moreover, A 2 = (^a w a) As a 
consequence, A has eigenvalues dbyj^, except in two dimensions, where 
A is proportional to the identity. Then, the determinant of M can be evalu- 
ated in the basis of eigenvectors of A, and in all cases, detM can be seen to 
vanish if A = ° . For A = 0, the determinant can only vanish if u>„ = 0. 

cosh 2 6» ' J a 

Otherwise, it can only happen for A 6 R+, which proves strong ellipticity 
in C — R+, and any even dimension. This completes the proof of Theorem 
O □ 
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5 Heat kernel in an infinite cylinder 



In what follows, we present the heat kernel for (P 2 ,B) in an infinite cylinder 
M = R + x j\f of any even dimension. By cylinder we mean that the metric 
is of the type ds 2 = dx 2 m + ds^, where dsjj is the metric of the closed 
boundary j\f. 

In order to determine the heat kernel, it is useful to note that the chiral 
bag boundary conditions in equations (|4.1j) and l|4.2|) are equivalent, for each 
eigenvalue of the tangential part B of the operator P, to Dirichlet boundary 
conditions on part of the fibre, and Robin (modified Neumann) on the rest. 

In fact, let's first notice that the operators V+ = U+ !\ + n and V- = U ~ U $~ 

' r T coshr 9 costr 9 

are self-adjoint projectors, and they satisfy V+ + V- = 1 splitting V into 
two complementary subspaces. 

Let £ = x m -x m , and 77 = x m +x m and as before, let y = (yi,?/2, — ,2/m-l) 
be the coordinates on the boundary and x = (y,x m ). If we call <&j(y) 
the eigenspinors of the operator B = j^mjada (with a = 1,2, ...,m — 1) 
corresponding to the eigenvalue u, normalized such that 

J2^(y)My') = s m - 1 (y-y) 

with S m ~ l the Dirac delta function, and 

J dy <j>u(y)<t>uj(y) = 1 j 

dM 

we can expand ip(y, x m ) = ^ w fu(x m )(f>u(y). If tp = V+ip, then the condition 
(|4.1jl is identically satisfied, and only (|4.«S|) must be imposed at the boundary 
which, for each oj, reduces to 

cosh 0e~^ (d m + oj tanh 0) f u = . 

Since the factor to the left of the parenthesis is invertible, this is nothing 
but a Robin boundary condition. 

In the subspace ip = V-ip, the boundary condition (|4.1|l reduces to 

coshfle^/o, = 0, 

while (|4.3|) requires 

uf u = 0. 

Thus, in this subspace, both boundary conditions are nothing but homoge- 
neous Dirichlet ones. 
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As a consequence, the complete heat kernel can be written as a Dirichlet 
heat kernel on V-V and a Robin heat kernel on V+V . For the convenience 
of the reader we make the single ingredients explicit j^2] an d write 



Kit; x, x) 



K(t;x,x')(V-+V^ 



1 



+ 



'4mt 
1 

/ 4TTt 



At 



-I -v 

e At + e « 



-1-2 v ,t/ ^ (aali .^"^^-^^^er/c^^,*)]} P 4 



1 



+ - 



Mvrt 

2ru m 



4t — g 4t 



cosh^(0) 

where u w (i],t) = — \/tu tanh(#) , and 



£>£(y)^(y)e-^'|(e 

[l + y^rt)o;tanh^ ( ' ? '* )2 er/c[^(7 ? ,t)] 



(5.1) 



g 4t 



2 Z" 00 
erfc{x) = —= I d£e 
V 71 " Jx 



is the complementary error function. Note that (J5.1|) is a direct generaliza- 
tion of the heat kernel given in [H] for the two-dimensional case, which, in 
turn, coincides with the Fourier transform of equation (101) in 4 for an 
antiper iodic boundary fiber. 



6 Meromorphic properties of the zeta function 

Let us now analyze the boundary contributions to the global zeta function 
related to (|5.1j) . We first note that global quantities are necessarily divergent 
due to the non-compact nature of our manifold M = R + x J\f . This is not a 
severe problem because the result 1)5. 1JI allows us to identify easily the bulk 
term leading to a divergent contribution when integrated. In particular, it 
is the first term in 1)5. lj) that represents the heat kernel on the manifold 
R x J\f . In the following, without changing the notation, we will first ignore 
this term and this will allow us to determine the boundary contributions to 
the global zeta function. Alternatively, as we will show afterwards, one can 
introduce a localizing function of compact support in ()5.1j) such that the 
trace gives a finite result. 
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Let us consider the trace of ()5.1|) ignoring the first term. It is convenient 
to perform the Dirac trace (tr) first. Since 



/ 2n + n* 



2 m , 



Sosh 2 (9), 

the trace of the 'boundary' heat kernel reduces to 



TrK(t;x,x) = — y^utanhfle" 



x / dx m er/cK(2x m ,t)]e^ m+ ^( 2:c - i ), 



where the second and third term in (|5.1j) have cancelled each other. Now, 
using that 

' -^A+^(2W) er/cM2w)] " 



1 d 



2 dx r 



^ + lo tanh 6 e u ^ 2Xm ^erfc [ Uuj ^x m , t)} 

Tit 



we get 



7 r 
TrK(t; x,x) = —~Y^ ^ ' e uM erfc K(0, t)] - 1 



Ve 



e cosh^ e 



1 + erf(ujVt tanh ( 



-art 



Here we used erf(x) = —erf(-x) = 1 — erfc(x). 

Now, we can Mellin transform this trace, to obtain the 'boundary' zeta 
function of the square of the Dirac operator in the infinite cylinder 



4T(s) 
2™ 



4T(s 



E 
E 



dt t 



s-l 



g cosh^ 



-U) Z t 



dt t s ~ 1 e'^eerf(uVi tanh 0) 



= Cl(s,P 2 ) + C2(s,P 2 ). 

The first contribution can be readily seen to be 

Ci( S ,P 2 ) = ^(cosh 2 ^-l)C(s, J B 2 ), 



(6.1) 



(6.2) 
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where B is the operator defined in Section [5J 

As for the second contribution to ()6.1|) . it is given by 



2™ „ f°° , -J>t 2 /-(^v^tanhe) 



After changing variables according to y = ^ c °r^f ; and interchanging 
integrals, one finally gets 

C 2 ( S)J P 2 ) = ^fe^ cosh 2 ^^>gnH (^P 



x — / ^(l + y 2 )" 5 " 



/7T 



^cosh 2s ^(2s, J B)— / c%(l + y 2 ) s 5 



sinh ( 







ll'(.s) ,v '^tt 

sinh6>cosh 2s 9rj(2s,B) 



i r( s + i) 



2^ r(s) 

x 2 Fi Q,^ + S ,|-sinh 2 ^ . (6.3) 

The structure of the zeta function is similar to the structure found for spec- 
tral boundary conditions, see e.g. [Hlj- In particular, the analysis of the 
zeta function on M has been reduced to the analysis of the zeta and eta 
function on M. 

As already commented, from ()6.2|) and (|6.3|) one can determine the po- 
sitions of the poles and corresponding residues for the zeta function in any 
cylindrical product manifold, in terms of the meromorphic structure of the 
zeta and eta functions of the boundary operator. For the rightmost poles, 
explicit results can be given in terms of the geometry of the boundary Af. 
For example, for s = {m — l)/2 we see that 

ResC^^-,^) = i(cosh-^-l)Resc(^, J B 2 ) 

= i (cosh- 1 0-1) I r ' 2"Vo«AQ. 

Because £2 does not contribute, given rj(2s,B) is regular at s = (m — l)/2 
|21j . this equals Res £(s, P) and is the result expected from the calculation 
on the ball [23! ■ ^ or $ = the residue disappears as is known to happen 
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for the standard local bag boundary conditions 21J. Further results can 
be obtained by using Theorem 4.4.1 of [U- Given we considered the case 
without 'potential', it is immediate that 

Res C (^,P 2 )=0. 

Also, for the particular case of s = 0, the fact that £(s, B 2 ) and r/(2s, B) are 
regular at s = shows that C(0,P 2 ) = C(0,P) = 0. 

Given the local heat kernel (|5.1|) . a local version of the results of this 
section is easily obtained. To this end, we use a localizing function with 
compact support near the boundary, such that its normal derivatives at the 
boundary vanish, 



f(y,x m ) 



dx 



0, n € IN. 



Furthermore, we let P 2 denote the operator P 2 on the double R x M of 
R+ x M, and we extend the localizing function as an even function to the 
double. We use the notation / = f(y,x m ), fj\f = f(y,x m = 0) and / 
for / on the double. Introducing the local versions ((f , s, P 2 ), ((f,s,P 2 ), 
C(fj\r, s, B 2 ) and r/(/_A/, 2s, B) of the zeta functions and the eta function, 
(|5.1j) and previous calculations show that the following theorem holds: 

Theorem 6.1 

r(s)c(f,s,p 2 ) = r(s)^af,s,P 2 ) 

A ( cos h 2s 9-l)((f M ,s,B 2 ) (6.4) 



1 r(s + 1/2) 2s a r I 1 1 , „ 3 „;„i,2. 



4 

+ 20F^ r'oo"' sinhgcosh " s 2jFi ( §' § + s > I' " sinh " J ^> 2s ' 5 ) 

+h(s) 

where h(s) is an entire function. 

This result parallels Theorem 2.1 in ^Uj f° r spectral boundary conditions. 

7 Conclusions 

In this article we have considered the influence of generalized bag boundary 
conditions on the heat kernel and the zeta function. In order to guarantee 
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certain structural properties we have first shown the strong ellipticity of the 
boundary conditions. Work by Seeley |241 125j then shows the standard heat 
kernel expansion holds and so the zeta function can have only simple poles 
at s = m/2, (m-l)/2, ...,1/2, and s = -(2/ + l)/2, I € IN. This is the main 
result of this paper. 

Based on the strong ellipticity one might envisage the determination of 
the leading heat kernel coefficients for generalized bag boundary conditions 
as they are needed for the calculation of effective actions in gauge theories in 
Euclidean bags [3] . Special case calculations can serve to restrict the general 
form that coefficients can have, cylindrical manifolds providing a valuable 
example. Here, for P = i'jj'Vj, we have expressed the heat kernel and the 
zeta function of the associated second order operator on M. = R + x J\f 
in terms of the boundary data on N . In fact, this result, under certain 
restrictions, can be straightforwardly generalized to P = i'yj'Vj — (ft- In order 
that a separation of variables as presented succeeds we need d Xm = and 
{"fm, 4>} = {7,0} = 0. If this is satisfied, equations Q6.2|) and (|6.3|) remain 
valid, once the operator B incorporates the potential, B = 77 m (7 a V a — 
icft). We have thus a particular case involving a potential and Riemannian 
curvature and various restrictions on heat kernel coefficients will follow. 

A Appendix: 7-matrices 

Let m = 2m be the dimension of a Riemannian manifold. We denote by 
Jj , j = l,...,m, the self-adjoint 7-matrices projected along some m-bein 
system. These are defined inductively by 




starting from the Pauli matrices 




The 7-matrices satisfy the Clifford anti-commutation formula 




(A.l) 
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In the main body of the paper we will simplify the notation and we will not 
indicate the dimension explicitly In addition, we set 



l m +i = 7 = (-*) 7i-7m, 
which is the generalization of '75' to arbitrary even dimension. 
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